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(Dated: August 6, 2003) 

The Hamiltonian dynamics of two-component spherically symmetric null dust is studied with 
regard to the quantum theory of gravitational collapse. The components — the ingoing and outgoing 
dusts — are assumed to interact only through gravitation. Different kinds of singularities, naked or 
"clothed", that can form during collapse processes are described. The general canonical formulation 
of the one-component null-dust dynamics by Bicak and Kuchai' is restricted to the spherically 
symmetric case and used to construct an action for the two components. The transformation from a 
metric variable to the quasilocal mass is shown to simplify the mathematics. The action is reduced by 
a choice of gauge and the corresponding true Hamiltonian is written down. Asymptotic coordinates 
and energy densities of dust shells are shown to form a complete set of Dirac observables. The 
action of the asymptotic time translation on the observables is defined but it has been calculated 
explicitly only in the case of one-component dust (Vaidya metric). 

PACS numbers: 04.20.Dw, 04.20.Fy, 04.60.Ds, 04.60. Kz 



I. INTRODUCTION 

The theory of gravitational collapse and black holes 
has acquired many facets, extending from detailed astro- 
physical models of activity in galactic nuclei to the so- 
phisticated calculations of the black-hole entropy. Still, 
some long-standing central questions persist: Is the out- 
come of collapse always a black hole or can in some cases 
naked singularities be formed? Even if the cosmic cen- 
sorship conjecture will sometime be formulated as a the- 
orem, the final singularities inside black holes will require 
new physics. Hopes in curing singularities by a suitable 
form of quantum gravity are older than the concept of a 
black hole itself. Despite new advances towards a gen- 
eral formulation of quantum gravity, simplified models 
play still important role, not least because within such 
models one may address oneself more easily to a concrete 
physical question. 

The spherical gravitational collapse of null dust be- 
longs to this type of problems. Since Vaidya's discovery 
of an exact solution of Einstein's equations describing a 
"radiating spherical star" T] , null dust - also referred to 
in literature as "pure radiation field" - has been widely 
used as a simple matter source. Its energy-momentum 
tensor 



where kn is null, 



Taf3 = pkakp, 



(1.1) 



(1.2) 



may be interpreted as an incoherent superposition of 
waves with random phases and polarizations moving in a 



single direction. Null dust exhibits all features of the ge- 
ometrical optics limit of Maxwell's theory except for the 
polarization properties 0] , . In this limit the mass dis- 
tribution p can be identified with the square of the scalar 
wave amplitude. Alternatively, the energy-momentum 
tensor l|l.l|l can be considered as representing photons 
that move along null rays with the fiux vector determined 
by the amplitude and the null vector fc" 0. The light- 
like particles may be massless scalar particles or neutri- 
nos. The same energy-momentum tensor describes high- 
frequency gravitational waves. One can also reinterpret 
some exact solutions of Einstein's equations with null 
dust as spacetimes produced by zero-rest-mass fields. All 
these roles of null dust indicate that it is more closely re- 
lated to fundamental matter fields than ordinary dust 
from phenomenological massive particles. For more de- 
tails and references on the null dust and geometrical op- 
tics, see Appendix A in for the list of various exact 
solutions with null dust, see Q, and Appendix B in 
where some applications are also mentioned. 

Vaidya metric is an interesting model in respect of 
studying the singularities and their possible removal by 
quantum phenomena. Indeed, classical solutions with 
null dust contain both naked and "clothed" singularities. 
In the next Section 2 we shall briefiy review how spherical 
null dust solutions can be used to study the formation of 
naked singularities, or to describe mass inflation inside 
black holes. 

The main purpose of the present work is to construct 
the canonical theory of spherically symmetric spacetimes 
with null dusts. The general Lagrangian and Hamilto- 
nian framework which includes (single-component) null 
dust as a source into canonical gravity was developed in 
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Ref. 13- In Sec. 3, we generalize the Hamiltonian action 
given in |^ to include both radially ingoing and outgoing 
null dusts. We assume that the null dust components in- 
teract only gravitationally, i.e., even in the region where 
the two streams interpenetrate, the energy-momentum 
tensors of both components are conserved separately. In 
the same section, the action is reduced by spherical sym- 
metry. Since in this paper we are interested in the inter- 
action of outgoing and ingoing null dusts, let us add that 
recently the similarity solutions representing interacting 
null dusts (not necessarily with the same type of inter- 
action between the components as we have chosen) were 
studied ; and the exact static solution for two colliding 
spherically symmetric null dust streams in equilibrium 
was presented 0. 

Although the present paper does not address any quan- 
tum problem directly, our main motivation for this gen- 
eralization originates in some problems of the quantum 
theory of gravitational collapse. In Refs. and 
quantum gravitational collapse of a spherically symmet- 
ric null-dust thin shell has been studied (this is a kind 
of limit of the Vaidya metric) . The quantum mechanics 
constructed there is unitary and describes a bounce and 
re-expansion of the shell in spite of the fact that the cor- 
responding classical solution contains a horizon and a sin- 
gularity inside it. The problem mentioned above is: what 
is the nature of the quantum spacetime geometry around 
the quantum shell? The metric itself is not a gauge in- 
variant quantity in the sense that different gauges lead to 
unitarily inequivalent metric operators (see, e.g., 0). We 
have proposed, therefore, to study the quantum theory of 
a system consisting of two shells so that the second shell 
can probe the geometry around the first one and vice 
versa, carrying some information about it to the infinity 
where it may be deciphered by the asymptotic observes. 
The investigations of such a system JJiJ has revealed a 
rich set of states in which the shells cross. Such crossing 
can be described by a finite-step canonical transforma- 
tion in spite of the crossing process being instantaneous. 
The reason, of course, is that the shells are infinitesimally 
thin but carry finite energy and momenta. The canoni- 
cal transformation describing a crossing is rather involved 
and there is no obvious way of how it can be "quantized." 
(Usually, we quantize an infinitesimal canonical transfor- 
mation by making the function that generates it to an 
operator.) As yet, we have not managed in a direct way 
to find a suitable quantum version (i.e., a unitary op- 
erator) of the canonical transformation for the crossing. 
One hope, however, is that it could be found indirectly, if 
we study some limit of the system of two null-dust thick 
shells crossing each other: there must be a well-defined 
Hamiltonian for this process. Such a study is initiated 
by the present work. 

The gravitational part of the constraints can be sim- 
plified by Kuchaf 's canonical transformations introduced 
originally in the vacuum Schwarzschild geometrodynam- 
ics [Oi l and used subsequently for a null-dust thin shell 
model in (T3|. In Kuchaf 's procedure the quasilocal mass 



plays an important role. In case of colliding spherical lay- 
ers of null dust the quasilocal mass is well identified in 
Bondi-type coordinates used by Bardeen to study the 
effects of back reaction of the Hawking radiation. These 
coordinates will be used extensively throughout the pa- 
per. The Kuchaf transformation is carried out in Sec. 
4, based on a careful analysis of the resulting boundary 
terms. 

In Sec. 5, we reduce the action by a gauge choice so 
that the Hamiltonian for cross-streaming null dusts can 
be found. The gauge is shown to be regular, and the 
Hamiltonian is calculated. The result, Eq. H5.12|l . is non- 
local — one must be careful in calculating its variation. 
Then the dynamical equations implied by Hamiltonian 
(|5.12l) are obtained and, in Sec. 6, compared with the 
Einstein equations in the form that has been given by 
Bardeen [l3l |. 

In Sec. 7, we analyze the dynamical equations looking 
for integrals of motion. Two explicit integrals of motion 
per point are found. From them a complete set of Dirac 
observables is constructed. The transformation between 
the original phase space variables and the Dirac observ- 
ables so constructed is, however, not known in an ex- 
plicit form. Similarly, the expression of the Hamiltonian 
in terms of the Dirac observables is not known explicitly. 
Nevertheless, the transformation properties, as well as 
the physical meaning of the Dirac observables are avail- 
able. And the lack of the explicit algebra and dynamics 
of our Dirac observables does not obstruct our project 
of finding the quantum version of the shell crossing. We 
can work with the canonical coordinates of Sec. 5 and 
with the Hamiltonian H5.12|) . This will be done in the 
subsequent work. 

It turns out that the explicit transformation between 
the Dirac observables and the variables of the reduced 
theory, as well as the expression of the Hamiltonian in 
terms of the Dirac observables, can be found in the spe- 
cial case of one-component dust. This is shown in Sec. 
8. To achieve this aim, we start from the extended phase 
space and constraints (before the reduction by a gauge 
choice of Sec. 5). We use the Vaidya coordinates as 
our covariant gauge condition and introduce the corre- 
sponding embedding variables. The coordinates in the 
extended phase space are transformed to the embedding 
variables and the Dirac observables. Then the reduction 
to the physical phase space spanned by the Dirac observ- 
ables is carried out and all desired transformations are 
written down explicitly. In this_way the action is trans- 
formed into the Kuchaf form |j9| — cf. Eq. (|H?TH1) . This 
section can also illustrate how a nontrivial dynamics of 
Dirac observables based on the asymptotic symmetries 
can be constructed. These ideas are explained thoroughly 
in Ref. 

The canonical theory of one-component spherical null 
dust, starting from the general formalism of Q and 
Kuchaf 's transformations, was recently developed by Vaz 
et al |l5j |. The authors also discussed the quantization 
procedure. However, their final form of the action is not 
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in the Kuchaf form and their theory is not gauge invari- 
ant. 



II. NULL DUST AND FORMATION OF NAKED 
SINGULARITIES 



As we have aheady stated in the Introduction, the null 
dust models have been used to clarify, in classical con- 
text, the formation of naked singularities during a spher- 
ical gravitational collapse. It has been known for about 
30 years that during spherical collapse of ordinary dust, 
the shells of dust may cross under suitable condition to 
generate the density singularity outside the event hori- 
zon (see e.g. for a review). Since it is possible to 
determine the motion of dust through the singularity, 
this "shell-crossing singularity" is considered as relatively 
innocent. More serious is the "shell- focusing singular- 
ity" which may form at the center of a collapsing spheri- 
cally symmetric inhomogeneous dust cloud described by 
a Tolman-Bondi metric. When collapse to the center pro- 
ceeds sufficiently slowly, a "past-null" singularity which 
is visible from infinity may arise at the center jl7l | . Such 
singularity is a "strong curvature singularity" , as defined 
e.g. in [ig. 

That a shell focusing naked singularity may form in 
case of imploding null dust, was first noticed by Hiscock 
et al , and analyzed in detail by Kuroda |23| and Pa- 
papetrou [2l| . HoUier proved that this singularity is 
a strong curvature singularity. A number of more recent 
papers analyzed this problem within classical theory . 
To describe the results and to introduce in more detail 
the following sections, we first need a few facts about the 
Vaidya metric. 

The energy-momentum tensor is the "primary" 

measurable quantity. Since the lightlikc normalization 
(|1.2|) is preserved by an arbitrary scaling of ka , 



Schwarzschild curvature coordinate R: 



k = Afc", A > 0, 
we can simultaneously rescale the "density" p, 



(2.1) 



(2.2) 



so that the form (|1.1|) of the energy-momentum tensor 
is preserved. By taking A = p^^^ we can eliminate the 
scalar p and give the energy-momentum tensor l|l.l|l just 
in terms of a single null vector, 



in the form 



(2.3) 



(2.4) 



This choice maximally simplifies the canonical descrip- 
tion of null dust. It will be employed in the following. 

The simplest form of the Vaidya metric describing 
spherical implosion of null dust is achieved by us- 
ing the incoming null coordinate V, together with the 



ds^ = - 1- 



2M{V)\ 



dV^ + 2dVdR + R^dif. (2.5) 



Here M{V) - the "advanced mass", measured at past 
null infinity T~ - may be an arbitrary, non-decreasing 
function. In the time-reversed case of outgoing spherical 
cloud of null dust, the simplest description is provided in 
terms of the outgoing null coordinate U: 



ds' 



2M{U) \ 
R ) 



dU^ - 2dUdR + R^dVt^, (2.6) 



where the "retarded mass" M(U), measured at future 
null infinity X"*", may be an arbitrary, non-increasing 
function. For M = constant, the metrics (|2.5ll and 
()2.6|l go over into the Schwarzschild metric in the ingo- 
ing, respectively outgoing Eddington-Finkelstein coordi- 
nates The metrics remain simple if, instead of null 
coordinate V, respectively U, Eddington time-coordinate 
= V — R, respectively = U + R, is introduced. 
However, Vaidya's metric becomes more complicated in 
the Schwarzschild-type coordinates 

In case of ingoing null dust the metric H2.5|l together 
with the Einstein equations imply 



where My = dM/dV > 0, 



Analogously, for outgoing null dust we get 



"'3 47ri?2 
where M^u = dM/dU < 0, 



h.+ t.+ — 1+1+ 

'^a '^0 ~ 'q '/3 I 



47ri?2 



1/2 



(2.7) 



(2., 



(2.9) 



(2.10) 



Consider a spherical layer of null dust imploding at 

V = Q from I~ (see Figs. 1 and 2). In region y < 
spacetime is flat. For > the gravitational field is 
given by Vaidya's metric l|2.5|l . the energy- momentum 
tensor by Eq. igSl). The mass function M{V) — for 

V < M{V) > is non-decreasing for V > Q. The 
scalar invariant quadratic in the Riemann tensor is equal 
to 48M2(T/) /R^ - a naked singularity may arise at i? = 0, 

V -+ 0-^-. If there exists an outgoing radial null geodesic 
which starts at R — 0, V — 0+ and ends at I+, the sin- 
gularity is globally naked. The analysis of such geodesies 
is not easy - Kruskal-type ("double-null") coordinates 
which would enable one a suitable description of both 
incoming and outgoing null geodesies can be introduced 
explicitly only when M{V) is a linear or exponential func- 
tion (2Jj. In the linear case, M{V) = aV,a = constant, 
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FIG. 1: The Penrose diagram of a thick spherical layer of null 
dust collapsing from T~ and forming a Schwarzschild black 
hole. The region containing null dust is shaded, the apparent 
horizon is indicated by dot-and-dashed line, the event horizon 
by full line. 




FIG. 2: The Penrose diagram of a thick spherical layer of null 
dust collapsing from X~ and forming a global naked singu- 
larity. The apparent horizon is indicated by dot-and-dashed 
line, the event horizon by full line, and the Cauchy horizon 
by dashed line. (Based on Fig. 3 in |20|.') 



it can be proven [20j that the singularity will be glob- 
ally naked if a < 1/16 and null dust inflow will stop at 
some V — Vq so that spacetime becomes Schwarzschild 
with mass M{Vq). This case is shown in Fig. 2. In 
the cases when constant parameter a > 1/16, the sin- 
gularity is hidden inside the event horizon - a black is 
formed as illustrated in Fig. 1. Therefore, the forma- 
tion of a naked singularity indeed requires a more gentle 
{a < 1/16) inflow of matter than the formation of a black 
hole {a > 1/16). Recently, naked singularities in higher 
dimensional Vaidya spacetimes were also considered . 

An interesting application of both ingoing and outgo- 
ing charged Vaidya metrics appears in the work of Pois- 
son and Israel |2Q| on internal structure of black holes. 
An initially static Reissner-Nordstrom black hole is per- 
turbed (non-linearly) by cross-flowing outgoing and in- 
falling null dusts. The region where the streams intersect 



lies below the outer horizon. We refer a reader in partic- 
ular to Fig. 3 in [2Q] in which - as in our Fig. 3 in the 
following - there are four regions with different values of 
the mass parameter outside cross-flowing streams. Pois- 
son and Israel show that the mass parameter inflates to 
classically arbitrarily large values on the Cauchy horizon 
inside the black hole. The timelike Reissner-Nordstrom 
singularity, which is locally naked, is thus probably 'con- 
verted' into a null or spacelike singularity. The mass 
inflation and the instability of the Cauchy horizon inside 
black holes have been also studied by other authors [27l |. 

The mass inflation represents another classical effect 
which, in a quantum context, requires a canonical theory 
for cross-streaming null dust flows. We now turn to its 
formulation. 



III. CANONICAL FORMALISM FOR 
SPHERICALLY SYMMETRIC SPACETIMES 
WITH TWO-COMPONENT NULL DUST 

We shall start from the geometrodynamical approach 
to a spherically symmetric geometry as described in |llj . 
Up to very few exceptions (noted in the following) we 
adopt the same notation and natural units G = c = 1. 

In geometrodynamics, a spherically symmetric space- 
time is generated by evolving a spherically symmetric 
three-dimensional Riemannian space S with the line el- 
ement described by two functions, A(r) and R{r), of a 
radial variable r: 

da^ ^ A2(r)dr2 + R^{r)dn''. (3.1) 

Here dVl^ = dO^ + sin'^ 6dip^, -R(^) is the curvature radius 
of the 2-sphere r — constant; under transformations of r, 
R{r) behaves as a scalar, A(r) as a scalar density. The 
range of the radial variable r depends on the problem 
we wish to discuss. If null dust falls into a primordial 
Schwarzschild black hole (which thus existed before an 
"irradiation" by dust), variable r g (— cxd,+oo) covers 
both ("right" and "left") infinities, whereas if an infalling 
"ball" of null dust creates a Schwarzschild black hole, the 
point r = is the center of the ball, and r G< 0, -foo). 
As discussed in the Introduction, our main interest will 
be in treating the thick shell of null dust falling in or out 
with Minkowski space inside. Then again r £< 0,-|-oo). 
Clearly, appropriate boundary conditions must be for- 
mulated and boundary terms analyzed depending on the 
physical situation considered. 

If a spherically symmetric spacetime is foliated by 
spherically symmetric leaves E, labelled by a time pa- 
rameter t G R, the metric functions A and R depend 
also on t. The leaves are related by the lapse function 
N{r, t), and by the only non- vanishing component of the 
shift vector N^{r,t). Modulo boundary terms, the vac- 
uum dynamics of the spherically symmetric gravitational 
field follows from the canonical form of the ADM action, 
integrated over and If. It reads (cf. [ill, Eqs. (31)- 
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(33)): 



5*^ = J dtj dr{Pf,A + PrR - NH"- - iViJ^^), (3.2) 
where the gravitational super-Hamiltonian 



H'^ = R-^PrPa + -R-^KPl + A-^RR" 



(3.3) 

K-^RR'K' + ^A-ii?'2 - ^A, 



and supermomentum 

H'^ = PrR' - AP;, 



the dots and primes denoting the derivatives with respect 
to variables t and r. The momenta, Pr and Pa, obtained 
by differentiating the ADM spacetime action with respect 
to the velocities A and R, read 



Pa = ~N-^R{R^ R'N'' 



(3.5) 



Pr = -iv-i 

Inversely, 



A(P - P'TV) + P(A - (A7V)') . (3.6) 



A = ~NR-^{RPr- APA) + iAN''y, (3.7) 



R = -NR-^Pa + R'N\ (3.8) 

A general spacetime action for a one-component null 
dust has the form 3] 



=JdtJ d3a;L^^, 



where the Lagrangian is given by 

1 



L 



ND 



8tt 



with 



(3.9) 



(3.10) 



(3.11) 



i = 1, 2, 3 is a summation index. We inserted the factor 
{4tt)^^ into (|3.1U|) as compared with the general action 
introduced in U], so as to obtain simplest expressions af- 
ter the integration over 9 and ip is made. The symbol 
I" denotes the unique vector field defined by the relation 
Tfj,i, = Ifilv The fields represent maps from the space- 
time M. to the material space Z of the dust, each triple 
(Z^,Z^,Z^) being a name of a particular dust volume 
element. 

A general coordinate transformation in Z can be writ- 
ten as 

Z^^Z\Z\Z\Z^), det(^\^Q. (3.12) 



For Wi, we have Wi = (dZ^ /dZ^)Wj, so it transforms 
as a covector under the transformation H3.12|l . The La- 
grangian is manifestly invariant. Observe that the trans- 
formation group is Diff2^, an infinite-dimensional group, 
but it cannot be considered as a gauge group. The trans- 
formations H3.12|l are "rigid" or "global" in the sense that 
they are independent of spacetime points. 

In Ref . 3] , the transition to the Hamiltonian action for 
the null dust is performed with the result 



S 



ND 



dt / d-'x{PkZ''-NH'^''-N''Hi^'') , (3.13) 



(3.4) where 



H 



ND 



gabHNDjjND 



and 



ND 



^PkZ' 



(3.14) 



(3.15) 



the variable Pk is the canonically conjugate momentum 
to Z'^ and it transforms as a covector with respect to 
(Em; g'^^ is the inverse of the metric gab induced on the 
leaves S by the spacetime metric. 

During the transition to Eqs. (|3.13(l - (|3.15(l . important 
relations have been obtained and used. Some of them 
cannot be derived from the reduced action H3.13|l any 
more. Thus, it has to be supplied by the formula for 
the vector field P which follows from the Hamiltonian 
equations (see Ref. 3], Eq. (6.29)-(6.31); notice that on 
the right-hand side of Eq. (6.30) the sign should be "-f-"): 



1° 



9 



ah 



-.Y° 



where n" is the unit future-oriented normal vector to 
the foliation and Y" is determined by the description 
Y" = Y" {x^ , x'^ , x^) of the foliation in the coordinates 
Y" of the spacetime and x^ of the foliation surfaces. The 
abbreviations 



W 



g^^ =g^'z:^Z^,, \g\=detgab 



are used; g^^ is the induced metric on E expressed in the 
basis Z\. 

To reduce action H3.13|l by spherical symmetry, we 
choose 

{Z\Z',Z') = {^,9,^) , 

where $ labels the comoving spherical surfaces of the 
dust while 6 and ip mark the radial lines of the dust. From 
the symmetry, we have 9 = ip — 9' — ip' — 0; we assume 
that the spacetime coordinates 9 and ip are identical with 
the dust coordinates. The conjugate momenta are (cf. 
Eq. (6.18) of 0) 

(Pi,P2,P3)EE (n,o,o) . 
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For the metrics, we have 



f 2A^N''dtdr 
A^dr^ + R^dn^ 



(3.16) 



and 



Then, 



(7Vi,7V2,A^3) = (^'',0,0) . 
g,j = diag(A2$'-^ R^, sm^ d) 

|n$'| 



w = 



A 



and 



r 



AR 



(3.17) 



where e = sgn(n(f>') describes the two possibilities of ra- 
diaUy out- and ingoing dust. Thus, we obtain after inte- 
gration over the angles: 



= I dtdr 



n|*-(f 



(3.18) 



As 500 < 0, we must have N / A > N^. If $ is chosen to be 
positive, we must have e$' > 0. For ingoing dust, $' > 0, 
while for the outgoing one, $' < 0. Hence, e = +1 for 
the in- and e — —1 for the outgoing dust. 

The generalization from one to several dust compo- 
nents that do not directly interact with each other (indi- 
rectly means through the gravitation) is simple: several 
terms in the action similar to the one-component action 
must be added. This, for two components, leads to the 
null-dust action of the form: 



:<ND 



dt dr 



7V(n+$' 



f |n_$'_ 
-n_$' ) 



The total action is 



where 
H 



dtdr{PAA+PRR+Il+^++n^^^-NH-N''Hr) , 

(3.19) 



PrPa 
R 

|n+<i>; 



^ _ 
|n_$'_l 



RR" RR'A' 



A2 



R^ 
2A 



A 
2 



A 



AP; 



(3.20) 
(3.21) 



A 

PrR' a I ^^+^+ 

Eq. H3.17|l must be written for each dust component 
defining in our case the vector fields and P, so that 



T 



af3 



47r V + 



(3.22) 



the total stress-energy tensor is the sum of the component 
terms because the total null-dust action is the sum of two 
terms of the form H3.10|l . 



IV. KUCHAR TRANSFORMATION 

The gravitational part of constraint H3.20|l can be sim- 
plified by a canonical transformation that replaces A by 
the quasilocal mass M as it has been first shown by [ll| ■ 
One way to define the mass is to use the coordinates 
introduced by Bardeen The metric reads 

ds^ = -Fe^'^dV'^ + 2e'^dVdR + R^dQ^ (4.1) 

in the advanced Bardeen coordinates (ABC), and 



ds^ 



-Fe^'I'dU^ - 2e'l'dUdR + R^dVt^ (4.2) 



in the retarded Bardeen coordinates (RBC). Here, we use 
the abbreviation 



F =1- 2M/R , 



(4.3) 



and the pairs of functions M{V, R),tp{V, R) or 
M {U, R) , (j){U, R) specify the geometry. The function M 
is a scalar field (it is the quasilocal mass of the spheri- 
cally symmetric spacetimes) while ^ and are potentials 
determined up to the transformation 

ip i—f ^ + lnv' , (f) i—i- (p + lnu' , 

because the null coordinates U and V are determined up 
to 

V^v(V), U^u{U), 

u and V being two arbitrary functions of one variable (the 
prime denoting the ordinary derivative of a function of 
one variable). We shall work with both ABC and RBC 
but derive our formulae explicitly only for ABC. 

Since we foliate the spacetime by spherically symmetric 
spacelike leaves S, labelling the leaves by a time param- 
eter t G R, and use a radial label r on the leaves, we can 
substitute the foliation 

V = V{t,r), R = R{t,r) (4.4) 

into the line element H4.1|l and find 

ds^ = - (Fe^^y^ _ 2e'^VR)dt^ 

+ 2{-Fe^'^VV' + e'''VR' + V'R)dtdr (4.5) 
+ {-Fe^'I'V'^ + 2e'f'V'R')dr^ + R^dVi^, 

where the dots and commas again denote d/dt and d/dr. 
Comparing H4.5|l with the general ADM form of the 
spherically symmetric line element (|3.16ll , we can express 
function A, and the lapse and shift functions in the form: 



A^ = -Fe^'^y'^ + 2e^V'R' 



N 



e^{VR' ~V'R) 
A ' 



(4.6) 



(4.7) 
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-Fe^HV + e-^VR' + e^V'R 
N = X5 • (4-8) 

Returning now back to the expression 1)3.51) for the mo- 
mentum Pa, we may substitute for A, N, and to ob- 
tain 

Pa = fli . (4.9) 

From the last equation we can express V regarding (|4.t)|) 
in the form 

-eV' = -^ + |. (4.10) 

From Eq. I|4.6|l . we can calculate function F, Eq. 1)4.3)1 . 
in terms of the canonical data: 

Substituting for F into Eq. ()4.1U)I . we obtain e^V also 
entirely in terms of the canonical data. Finally, combin- 
ing Eqs. ()4.3)l and 1)4.11)1 we can express the mass M as 
a function of the canonical data, 

M = -R^^Pa - -A-^RR'^ + -R. (4.12) 
2^2 2 

Notice that the equations 1)4.11)1 and ()4.12)) are of exactly 
the same form as in the vacuum Schwarzschild geometry, 
cf. [13. They can also be interpreted in a similar way. 

Clearly, the corresponding expressions in RBC are ob- 
tained by replacing 

everywhere, and by changing the sign of N (keeping N 
positive) . 

In vacuum spherically symmetric spacetimes one can 
find a canonical transformation to a new set of vari- 
ables in terms of which the constraints remarkably sim- 
plify p^ . The transformation replaces the canoni- 
cal variables {A,Pa,R,Pr) by new canonical variables 
(M, Pm, R, Vr) where M is given by Eq. 

Pm = R-'F-'APa, (4.13) 



i? = i?, (4.14) 



■Pr = Pr- ^R-'AP^ - -R-^F-^AP^ 

- R-'A-^F-^ [{APAYiRR') - {APA){RR'y] , 

(4.15) 



where F is given in terms of "old" canonical variables by 
()4.11l) . This transformation was shown to be canon- 
ical without employing constraints or dynamical equa- 
tions. We can make easily sure that it remains a canon- 
ical transformation also on the extended phase space 
which includes the null dust canonical variables $_|_,n_|_ 
and $_,n_. 

The only point we have to check is that the bound- 
ary terms vanish as in the case of primordial vacuum 
Schwarzschild black holes analyzed in 11]. We con- 
sider the following arrangement: An outgoing thick shell 
of dust starts at the past singularity and runs through 
the spacetime to reach X+ ; another thick shell starts at 
1^ , crosses the first one and reaches the future singu- 
larity (Fig. 3). Under both shells, the spacetime is fiat 
(Schwarzschild mass parameter M = 0). In the future of 
both shells, we have Schwarzschild spacetime with mass 
M = M3, in their common past M — Mi, and above 
both shells, M — M2. The hypersurface S is assumed to 
run from the regular origin i? = in Minkowski space- 
time to the i° of the Schwarzschild spacetime above the 
thick shells, as indicated in Fig. 3. 




FIG. 3: The Penrose diagram of two crossing thick spherical 
layers (shaded) of null dust. One layer is collapsing from 
I~ into a future Schwarzschild singularity of mass M3. The 
other is expanding from the past singularity of mass Mi into 
T"*". Below the layers the spacetime is flat (M = 0), above 
the layers the Schwarzschild mass is M2. A general Cauchy 
surface E going from the origin R — through the crossing 
region to is shown. The past and future event horizons are 
indicated by dashed lines. 

The ABC near i? = of S must be a transform of the 
advanced Minkowski coordinates Vm, R,0,ip'. 

Vm = viV) , 

where the function v could in principle be determined by 
boundary conditions at infinity and integration of Ein- 
stein's equations through the thick shells. For Vm, R, 6 
and if , the Bardeen potential ip is zero, hence 

V'(V) = Inv'iV) . 

At X~, we assume V to be the advanced time. Thus, 
V is the Eddington-Finkelstein coordinate above and in 
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the past of the thick shells. Within the ingoing shell and 
above the outgoing one, V is the advanced Vaidya coor- 
dinate. We can, therefore, set ip — in a neighbourhood 
of X^, and in the part of X"*" lying above the outgoing 
shell. 

We have thus to pay attention to the behaviour of the 
canonical variables A, R, Pa and Pr at the origin and 
at spatial infinity. Since we assume that thick shells of 
radiation do not extend to spatial infinities, there our 
fall-off conditions are identical to those considered in ITl| . 
Therefore, we assume the canonical variables at infinities 
to behave as follows: 



A{t,r) = 1 + M2it)\r\-' 
R{t,r)^\r\ + 0{\r\-'), 
PA(t,r) = 0(|r|-), 
Pfl(t,r) = 0(|rr(i+^)), 



Oi\r\ 



-(! + £ 



(4.16) 



where 0(r~*) means that this term falls off as r~*, 
whereas its j-th spatial derivatives as r"*^*"'"-'); the ab- 
solute values appear since in the case of primordial black 
holes r g R. The mass term can depend on t. The 
Schwarzschild time T outside the thick shells along a slice 
S is assumed to behave as 



(4.17) 



whereas the Lagrange multipliers in the action, the lapse 
and shift, satisfy conditions 

N ^N±{t)+0{\r\-'), ^0{\r\-'). (4.18) 

Since advanced time V is given in terms of canonical 
coordinates by Eqs. (|4.1U|) . (|4.11|) with = 0, the fall- 
off conditions 14.16|l imply 



V^\r\+2Moo\n 



(— 

V2Mo 



+ Too{t) + O{r-'), (4.19) 



so that 



V' = l + 2Moo-n + Oi\r\^^-'). (4.20) 



At the origin of Minkowski space the boundary conditions 
are the same as in 0. The 4- metric is flat, the 3- metric 
must be smooth. The hypersurface S must meet the 
origin parallelly to T = constant hypersurfaces to avoid 
conical singularities. Therefore, we require at r ^ 



T'(r)=0, R = r + Oir^), 



(4.21) 



evaluated at the boundaries, can also be written in terms 
of the advanced time variable V and function F as 



B = -RSRln 
2 



V'e'f'FR - 2RR' 



e^V'FR 



(4.24) 



where V'e'l' = V^^ &i R = Q and V' e'l' = V at R oo. 
Regarding the boundary conditions H4.16|l - H4.22|) . which 
imply 5R ~ 0(|r|~^) at infinity(ies) and SR ^ O(r^) at 
the origin, we easily see that the boundary terms vanish 
in all the situations we are dealing with. 

The constraints (jX^ and acquire now much 

simpler form: 

AH = -F-Hl'R'-FFMVR+\^+<S>',\ + \n^<i>'^\, (4.25) 



Hr = VrR' + PmM' + n+$V + n_$'_ , (4.26) 

where F is again given by (|4.3|l . but M{r) and R{r) are 
now understood as new canonical variables. (Notice that 
in (|4.25|l we scaled H by A.) From Eq. (|4.1U|) we find 
the advanced time coordinate in terms of new canonical 
variables also in a simplified form: 



"/-y p,, - p-^R'. 



(4.27) 



As we have seen above, we may assume that n_<i>'_ > 
for the ingoing and n+<i>^ < for the outgoing dust. We 
then combine Eqs. (|4.25|) and H4.26|l to obtain: 

Hr±AH^ {R' T FPm){Vr T P-'M') + | ^JJ^ Jr ■ 

(4.28) 

The total action (|3.19|) - (|3.21|) can now be written in the 
form 



S = 



dt dr 



PmM + VrR + n+$+ + n_$_ 



where 



N+H+ - N^H^ - 7V(oo)M(oo)] 



{R' + FPm){Vr + P-^M') + 2n+$', 



(4.29) 



+ ' 



= {R! - PPm){Vr~ P-^M') + 2Ii^<^'_ 



and 



Fm = r + P = To + r, V = 1. (4.22) 
The boundary terms, given by the expression (95) in [T]|. 



2 \ A 



2\ A 



(4.30) 
(4.31) 



B = -RSRln 
2 



RR' + APa 



RR' - APa 



(4.23) 



The action has also been supplied with the ADM bound- 
ary term. 
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V. REDUCTION BY A GAUGE CHOICE 

In this section, we shall reduce the action H4.29|l by 
choosing a gauge and by solving the constraints. The 
equations will simplify strongly, and this will help us to 
perform further transformations. 

Let us try the gauge defined by the conditions 



R = , Pi 



M 







(5.1) 



We have to check that this is a regular gauge, that is, 
the gauge functions defined by the left-hand sides of Eqs. 
(|5.1|) must have vanishing Poisson brackets with the gen- 
erators of gauge transformations. From the action H4.29|) . 
we obtain 



R^{N+~ N^)FPm + {N+ + N^)R' , 
Pm = [{N++N^)Pm]' + 



(5.2) 
(5.3) 



For the Poisson brackets it holds, of course, that 
{R, N+H+ + N-H-} = R and {Pm,N+H+ + N-H-} = 
Pm, where N+H+ + N-H- is a generator of gauge trans- 
formations for suitable values of and 7V_ . Regarding 
Eqs. H5.2|l and (|5.3() we see that these will be zero every- 
where at the gauge surface only if 



N++N_=0 



and 



F 



RF^ 



{N+-N-) = . 



These equations are equivalent to 



-i^Fexp {^2 J dR 
KF exp ( 2 / dR 



M' 
RF 
M' 
RP 



(5.4) 
(5.5) 



where K is a. constant. 

The constraints can be solved for Vr and M': 

M' = F(-n+$+ + n_$'_) . 



(5.6) 
(5.7) 



Eq. lEZI shows that M' > for F > while M' < 
for F < 0. Moreover, the integrands in Eqs. (|5.4|l 
and (|5.5|l are always regular and non negative and they 
converge for our arrangement: at i? = and near R — oo, 
there is no dust. Hence, both and N- will be non 
zero at i? = oo, and so the generators for which the 
Poisson brackets of the gauge functions vanish are not 
gauge generators. Hence, the gauge is regular. Eqs. (|5.4|) 



and (|5.5(l determine also the values of A''-!- and 7V_ for 
our gauge. Eqs. H4.30|l and (|4.31|) then imply 



TV 
where 



N 



h{R) 



2KAFexp U j\x ^ 



-n+<i>V 



(5.8) 



It follows that the gauge H5.1|l is i?-orthogonal one: the 
t-lines are the curves R — constant and the foliation is 
orthogonal to them. The foliation is determined by the 
second equation of H5.1|l , but the time function that is to 
be constant along the folios is not. We can choose the 
function, which will be called T, by imposing the condi- 
tion on N: Nipo) — 1. This determines the constant K 
so that the definitive N is 



N{R) = AFexp ^-2 j dx 



h{x) 



(5.9) 



Eq. H5.6|l together with the Hamiltonian equation for M 
and the above relations for the lapse and shift yield: 



FN 

M = — (n+$'+ + H_$'_) 



(5.10) 



Eq. H5.7|l can be considered as a differential equation 
for AI if the function F is written out according to Eq. 
(|4.3|) . This equation can be easily integrated: 



M{R) 



dxh(x) exp —2 



y 



(5.11) 



The Hamiltonian of the reduced theory is the value of 
the boundary term in the action H4.29|l that can be now 
calculated from the above equation with the result: 



dRh{R)exp [~2 dx 
\ Jr 



h{x) 



The reduced action is, therefore, given by 



(5.12) 



(5.13) 



The Hamiltonian H5.12|l is non-local. We can obtain 
the canonical equations for the Hamiltonian -ffred fol" 
lows. First, we vary iJfed with respect to h{R): 



SH, 



red 



dRSh{R)exp 



h{x) 



'2 / dx 
Jr X 

-2 dRu'{R)v{R) 
Jo 

where we have introduced the abbreviations 

h{y)\ 



uiR) 



J dx h{x) exp ^—2 J 



dy 



y 
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and 



v(R) 



r 

Jr 



dx 



5h{x) 



It is easy to express u{R) in terms of N{R) and M{R) 
using relations 15.11|l and (|5.9() : 



u{R) 



N{R)M{R) 

If 



We can see immediately that u(0) = and w(oo) = so 
that the last integral in 5i?i-ed is 



POO j*OQ 

-2 / dRu'{R)v[R) = - / dR 
Jo Jq 

Collecting all terms, we obtain 

^H,.cd = / dRSh{R) 
Jo 



2N{R)M{R) 5h{R) 



KF 



NjR) 
A 



R 



(5.14) 



The calculation of the canonical equations is now simple 
if we substitute for h{R) Eq. H5.8|l : 





= , 


(5.15) 




= , 


(5.16) 


X"-)' 


= , 


(5.17) 


^0' 


= . 


(5.18) 



n_ 
ri4 



The canonical equations (|5.15l) (|5.18|l together with 
Eqs. (|5.9() and (|5.11l) determine the dynamics of the dust 
and the metric. For example, Eq. (|5.10|) can be derived 
from Eq. H5.11|l by expressing M{R) first as a linear func- 
tional of ft,, similarly to the calculation of 5M in terms 
of 5h, and then by transforming h with the help of Eqs. 

The spacetime metric in the coordinates T, R, d and 
(p has the form, due to our choice of gauge. 



ds-" = -N^dT^ + K^dR^ + R^dQ. 



(5.19) 



The function A can be obtained from Eq. H4.(j|) . Indeed, 
the gauge condition R = r and Eq. (|4.8|l with A^*" — 
imply 

-Fe'^V + 1 = 0; 
this substituted into Eq. (|4.t)|) leads to 

1 



A = 



(5.20) 



where M{R) is given by Eq. (|OT|) . Eqs. (|^ and (|^?^ 

determine the function N . 



VI. REDERIVATION OF BARDEEN'S 
EQUATIONS 

The variation principle (|3.19|l can be checked by deriv- 
ing Einstein's equations in Bardeen's coordinates from it. 
The equations to be derived have the form [T3|: 



dtp 

dR = ^""^^^^ ' 
-dR - "^"^ ' 



(6.1) 
(6.2) 
(6.3) 



They are written in the ABC; analogous formulae hold 
in the RBC. 

To derive Eqs. (|6.1|l - (|6.3fl . we have first to express 
the stress-energy tensor using Eqs. (|3.22|l and H3.17|l . 
Transforming the metrics H4.1|) and H4.2|l into the coordi- 
nates T and R by the embedding relations V = V{T^ R) 
and U = U{T,R), and using the i?-orthogonality of our 
gauge, we obtain 



V = 



N 



eVVF 



U = 



N 



e't'VF 



and 



V 



-jf 



u' 



(6.4) 



(6.5) 



where the formulae for the RBC are also written out for 
later use; F is to be calculated from Eqs. H4.3() and l|5.11|l . 
The i?-orthogonality of the gauge implies, in the ABC, 

Y%^{V',l), n" = z.(l,0), 

where is a normalization factor. Eqs. (|4.1f) and H6.5f) 

yield 



1 



3'f'VF 







(6.6) 



Substituting this into Eqs. 13.22|l and H3.17|l . we obtain 
the right-hand sides of - 



Tbr — — 



1 



rpV _ 

" 2^ 

e^F^ 



F 

n+<i>' , 



(n+<i>+ + n_$'_ 



^ 4^i?2 

On the other hand, we also find the useful formula 



(6.7) 
(6.8) 
(6.9) 



t: 



^"^5 47ri?2 
and, in an analogous way. 

Hi $' 



T, 



Q/3 



47ri?2 



F^e'^'f'VaVf:! 
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Next, the derivatives with respect to ABC that oc- 
cur on the left-hand sides of Eqs. (|6.1|) - (|6.3|l are to be 
expressed in terms of the primes and dots. Clearly, we 
have for any function X{V, R): 

so that Eqs. (|6.4|) and H6.5() imply 



dV J N 



dRJy NVF 



(6.10) 

Analogous formulae for the RBC and X{U, R) read 



dU J R N 



dR 



Ny/F 



(6.11) 

Now, Eqs. 123, (ESI, (E3 and yield Eqs. 

(in3 and (ins. 

To derive Eq. H6.1|l . we start from Eqs. (|5.9() . H5.20|l . 
H6.4|l and obtain 



,^V = exp(^^2 p dx 



(6.12) 



Taking primes of both sides and using Eq. (|6.12|) again 
gives 

e^V^I^' + eH' = e'^V ^^^^^ 



R 



To calculate V', we substitute for V' from Eq. H6.5|l : 

\e^F) e'PF^ \ dV dV 
Of course, from Eq. (|4.3|l we have 
dF 2 dM 



dV RdV ' 

so that Eq. (|6.10() . together with the last three equations, 
entails 



a?/' _ 2 / 1 dM 
'dR" r[ ~ e'PF^'dV 



Then, Eqs. (EHI), W^^ - and (ES yield Eq. (^1)) . 



VII. DIRAC OBSERVABLES 

The equations of motion imply conservation laws, two 
for each dust component. We can use these conserved 
quantities as Dirac observables. 

Consider Eq. (|5.15|) : with the help of Eqs. H5.20|l and 
(El, it can be rewritten in the form 



1. _ e^FV^ ] - 
dT OR 



However, Eqs. H5.20|l . (|6.4|l and (|6.5|) also imply that 



Hence, $_ is conserved along T/- curves: 
'0^-\ 



dR 



(7.1) 



The same result, of course, follows from applying Eq. 

The second integral of motion is connected with the 
conserved current of Eq. (4.10) of It is interesting 
that it remains conserved even under the conditions of 
cross-streaming dusts. It has the form 



die'^Fn^] 
dR 







(7.2) 



To prove the equation, we express the left-hand side as 
follows: 



9(6"^ Fn_) 
dR 



n +e^F(^ 
OR Jv \dR 



and calculate the second term by means of Eqs. Hfi.lO|l . 
(1^:^ . and (|CT7|) : 



-n_ 



{e'f'FVy 

dR Jy'' e'f'Fv''''^ ^ e'f'FV 

We also obtain from Eqs. (|6.10|l and H6.4|l that 

'die''' FY 



n_ 



dR 



Then, Eq. H7.2|l follows immediately. 

Analogous conservation laws, this time along [/-lines, 
hold for the outgoing dust component: 



d^ 

dR 



1^™) (.3, 



These conservation laws can be utilized for the con- 
struction of a nice set of Dirac observables for the two- 
component dust. Let us define two functions C-(^) and 
f^-iV) by 

C_(l^) = $-, fi-{V) = e^Fn^ , 

as well as two functions C+{U) and ii^(U) by 

(+([/) = $+, ^l+{U)^e'^FU+ . 

The original null-dust canonical variables can then be 
expressed by means of these functions as follows 



<i>^{R) = C-{V{R)) , 
Il-{R) = ^l-{V{R))V'{R) , 

<^+{R)^C+iUiR)) , 
n+(i?) = -^,+ {U{R))U'{R) , 



(7.4) 
(7.5) 
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where we have substituted V or -U' for e^F and e'>'F 
with the help of Eq. 

Let us finally define the first pair w(C-) and w(C+) of 
the desired integrals of motion to be the functions inverse 
to C-(^^) and C+{U): 

C^{v{x))=x, (:+{u{x))=x. (7.6) 

The other pair and /i+(C+) is then defined by 

A*-(C-) = l^-{v{C-)) , M+(C+) = M+(«(C+)) ■ (7.7) 

The two pairs form a complete set of Dirac observables. 
We assume that the material space Z with coordinates C- 
and C.+ consists of two components, each with topology 

TIjl. 



and we obtain: 



C+G 



We also assume that the functions f(C-) and m(C+) are 
non-decreasing . 

Relations H7.4() and (|7.5() determine the transformation 
between (n„,<i>_) and (/i_,w) only implicitly, and so it 
is for the "+" variables. The function V{T, R) that oc- 
curs in them is, in fact, a complicated functional of the 
canonical variables. The reason is that the transforma- 
tion between the coordinates {V^ R) and (T, R) is solution 
dependent, and it can, therefore, be written as 

V = V[Il. (T, i?) , (T, R) , n+ (T, R) , (T, R) ;T,R], 

or 

y = i/[m-(C-), ^^(C-), M+(C+), "(C+); i?] , 

where the functions n_(T,i?), $_(r,i?), n+(T,i?) and 
$+(1", i?) in the first formula are considered as the initial 
datum for a solution (including the metric) at the Cauchy 
surface corresponding to a time value T. 

Formulae H7.4() and (|7.5|l can, nevertheless, serve to 
find the transformation properties and physical meaning 
of the Dirac observables. Indeed, Il± is a density in the 
material space as well as at the Cauchy surface T = con- 
stant. Divided by V or U', it becomes a scalar at the 
Cauchy surface and a density in the space with coordi- 
nates V or U. Thus 

/i+(c+)u(c+)dc+ , ^i-ic-m-)dc- 

are scalars with respect to diffeomorphisms in Z± and 
with respect to general transformations of the coordi- 
nates U and V. 

Let us express the total energy M in terms of the Dirac 
observables. Eqs. if^ . (|7^ and (|73|) yield 



M 



roc 

/ di?F(-n+$'|_ + n_$'_) 

Jo 



dR ^ + I dR 



U' 



V 



-4' , -V' ^C- 



With the correcting factors e~'^ or e""^, the integrands 
are just densities in the material space and the energy is 
a scalar. However, the explicit dependence of ^ or i/; on 
the Dirac observables is difficult to calculate. 

The physical meaning of the quantities ^i±{C±) can be 
made clear if we consider the arrangement of two crossing 
thick shells. Fig. 3. In a neighbourhood of X"*" that is 
intersected only by the outgoing thick shell, the solution 
can be written in the form of Vaidya metric that is given 
by Eq. ig^l) with = and M = M(t/). Then, Eq. if^ 
implies 



dM 

dc7 



U' 



Using Eqs. IjG.Sfl and H7.5|l with (f) ~ 0, we have 



dM du 

'du'dcl 



dM 
d^ 



and, similarly. 



dM 

dcI 



near for the ingoing thick shell. Hence, ^±(C±) are 
asymptotic energy densities in the material space. The 
corresponding variables w(C+) and w(C-) must then be 
defined by the retarded or advanced times U at X"*" and 
V at I^; u(C+) is the retarded time at which the mate- 
rial sphere C+ arrives at X"*" and similarly v(C-) is the 
advanced time at which the material sphere starts at 
I+. 

These simple results suggest that we may obtain more 
explicit formulae for the special case of one-component 
dust. This will be shown in the next section. 



VIII. VAIDYA SOLUTION 

Let us set H+ = so that there is only the ingoing 
component of null dust; from now on, we shall leave out 
the index "— " . The case with only the outgoing compo- 
nent can be treated in an analogous way. 

Eqs. l|HTTl and (|^ then imply that V' = tp{V) and 
so ip can be transformed away by a suitable choice of 
coordinate V. Moreover, Eqs. Ht).2|l and (|6.8|l show that 
M = M{V). Thus, we end up with the Vaidya metric. 

The canonical theory of the ingoing dust can start with 
the action (|4.29|) written out for H+ = 0: 



S = 



dt dr 



PmM + VrR + H$ 



N+H+ - N^H^ - 7V(oo)M(oo)] , (8.1) 
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where 

H+ = (R' + FPm){Vr + F-^M') , (8.2) 
H- = {R' - FPm){Vr- F-^M') + 2n<^' . (8.3) 

The variables in the action (|8.1|l can be transformed to 
embeddings and Dirac observables (Kuchaf decomposi- 
tion) as follows. 

Since Eqs. H4.6(l and 14.27|l together with ip ~ Q imply 



FP, 



M 



R' 



A2 
V 



(8.4) 



the first factor on the right-hand side of Eq. (|8.2|l is gen- 
erally non-zero. Hence the factor can be included in the 
Lagrange multiplier which we shall call N^: 



N+{R' + FP, 



M 



(8.5) 



The constraint equation = implies Vr + F ^M' = 
0; substituting this back into the constraint = 0, we 
get 



M 



IT 
F 



M' + n$' = 



(8.6) 



An equivalent action can thus be written as 



dt I dr PmM + VrR + n$ 
'- 

2iV_ {{Pm - R'F-^)M' + n$') 
N^iVR + P-^M') 



dtMoaTao , 



(8.7) 



where F is determined in terms of canonical variables 
by Eq. (|4.3|) . The Regge-Teitelboim boundary term 
in the action H8.ll) acquires the form MqoToo after the 
parametrization at infinity, as it is analyzed in detail in 
|llj |. (There would be another term of this type, corre- 
sponding to the left infinity, if null dust is collapsing into 
a primordial Schwarzschild black hole.) 

The action H8.7|l can be simplified by introducing new 
canonical momenta conjugate to M and R by 



n 



M 



Pi 



M 



P-^R' 



and 



TIr = Vr + F-HI'. 



(8.9) 



Substituting this into the Liouvihe form in H8.7|l . we find 



poo 

/ dr [Pm5M + Vr5R + n5$] = 

JQ 

/•oo 

= / dr[WM5M + J1r5R + J15^ 
Jo 

+ dr [F-\R'SM - M'SR)] . 
Jo 



(8.10) 



One can check by a short calculation that 

F-\R'SM - M'SR) = 

fRSRlnPy /R'RXnP 



.11) 



so that the transformation defined by Eqs. (|8.8(l and H8.9|l 
is indeed canonical, since with our boundary conditions 
the boundary terms resulting from the first term on the 
right-hand side of Eq. (|8.11|) drop out. 
Action H8.7|l then takes the form 



dt 



- 2A^_(nMM' + n$') - n^tvr 

dtMocToo- 



.12) 



The next step is to transform the variables M and IIm 
to V and fiy. Eqs. l|lT7jl with V = and imply 



Um = -V 



Let us define Ily by 



M(r) = - / dFny(?) , 
Jo 



(8.13) 



.14) 



SO that we have 



dr V' I d76nv (r) - MooST^, 
Jo 



= lim 



V{r) / dFmy(r) 





drV{r)5\lv{r)-Moo5To, 



Taking into account the boundary condition l|4.19|l and 
Eq. I|8.14|l . we have 



dr liv {r)5V{r) -5 dr V{r)I].v (r) 
Jo 



— lim 



2Moo In ■ 



2Af, 



6M{r) 



Since sufficiently large r always lie outside the matter, 
we can set M{r) = Moo within the limit so that the limit 
can be written as follows: 



lim 5 rMoo + 2 

r — ^oo \ 



dMMln- 



2M 
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Discarding all such total differentials in the Liouville 
form, we obtain finally 



Hence, 



drIlv{r)SV{r) 



Hence, the transformation (|8.13() and 1)8.14(1 leads to the 
action 



S= dt 



dr Hvl^ + HRi? + H$ 



(8.15) 



where — 2V'N-. The action has been transformed 
so that the original ADM coordinates A(r) and R{r) have 
been replaced by the embeddings V{r) and R{r). We 
observe that the ADM boundary term does not occur 
explicitly any more in the action H8.15|l : this seems to 
be a general property of embedding variables (another 
example is given in Ref. 0)- It seems that the role of 
the asymptotic time can be taken over by the embedding 
variables. The constraints in the action 1(8.15)1 are not 
identical with the conjugate momenta Hy(r) and Hfl(r) 
as it should be for an action after a Kuchaf decomposition 
0. The reason is, of course, that the remaining variables 
$(r) and H(r), which describe the physical degrees of 
freedom, are not Dirac observables and cannot have zero 
Poisson brackets with the constraints. In order to make 
the Kuchaf decomposition complete, we have to pass to 
the Dirac observables /i(C) and v{C). 

Eq. 1(7. 4() that relates H, $ and V to fi and v is valid 
for an arbitrary embedding V{r) and R(r): 

Hr) = CiVir)) , H(r) = /i(y(r))y'(r) , (8.16) 

but the variable V{r) plays now a different role in the 
formalism: it is itself a canonical coordinate in the action 
(|gTB|) . The problem that the function V{r) in Eq. ((7^ 
has not been explicitly known as a functional of canonical 
variables (of the reduced theory) has disappeared. 

Variations of the functions $(?"), ({v) and V{r) are 
related by Eq. ^J^ : 

dv 



dc ■ ■ n<i)' • 

H$dr = aCy'dr + li-^Vdr = ^xQdV + Vdr 

dv ' 



V 



The variations of the functions C(w) and v{C) are related 
by Eq. (|72J), 

5C + '^5v = , 
dv 



so that 



H$dr + Hyy = -jivdC, 
If we introduce the momentum Hy by 

H$' 



H$'\ 
Hy + — 1 Vdr 



Hy = Hy 



.17) 



then the transformation defined by Eqs. ((8.1t)() and ((8.17() 
is canonical. Substituting this into the action ((8.15(1 . we 
obtain finally 



S ■ 



dt 



(8.18) 



/ dr{nRR + UvV - N^Ur- N^Hv) 
Jo 



This action is in the Kuchai' form 
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